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Abstract 

The neutral Higgs sector of the minimal supersymmetric standard model (MSSM) in 
explicit CP violation scenario is investigated at the one-loop level. Within the context of 
the effective potential formalism, the masses of the neutral Higgs bosons are calculated at 
the one-loop level by taking into account the contributions of the following loops of ordinary 
particles and superpartners: top quark, the scalar top quarks, bottom quark, the scalar 
bottom quarks, tau lepton, the scalar tau leptons, W boson, the charged Higgs boson, the 
charginos, Z boson, the scalar and pseudoscalar Higgs bosons, and the neutralinos. Our 
calculation is an improvement in the sense that both the terms which are quartic in the 
electroweak coupling constants into account, and the pseudoscalar Higgs loop contribution 
are explicitly included. 
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1. Introduction 



For new physics beyond the standard model (SM) [1], the enlargement of the Higgs sector is 
considered as one of indispensable ingredients, especially for supersymmetric models. Among 
various supersymmetric extensions of the SM, the minimal supersymmetric standard model 
(MSSM) is the simplest one which has just two Higgs doublets in its Higgs sector [2]. In principle, 
in the models with multiple Higgs doublets, the violation of CP symmetry may be accomplished 
either spontaneously or explicitly by the mixing between the scalar and pseudoscalar Higgs parts 
[3, 4, 5]. In practice, however, it is known that at the tree level the Higgs potential of the MSSM 
conserves CP symmetry, because the complex phase can always be eliminated by rotating the 
Higgs fields. Thus, neither explicit nor spontaneous CP violation can happen in the MSSM at 
the tree level. Even at the one-loop level the scenario of spontaneous CP violation is excluded, 
because the radiatively corrected Higgs potcntcntial of the MSSM leads to a very light Higgs 
boson which is unacceptable by the CERN e"'"e~ LEP2 data. 

The remaining possibility at the one-loop level for the MSSM is then the explicit CP violation 
scenario. A number of investigations have been devoted to examine the explicit CP violation in 
the MSSM at one-loop level [6]. In those investigations it is observed that the tree-level situation 
is considerably modified when radiative contributions are included. In particular, the one-loop 
MSSM Higgs potential with explicit CP violation is found to allow the lightest neutral Higgs 
boson to possess a mass above the expemental lower bound from the LEP2 data. Therefore, 
the general consensus is that the MSSM at the one-loop level can accommodate explicit CP 
violation. 

Recently, Ibrahim and Nath [7] have investigated explicit CP violation scenario, paying their 
attention to the phenomenological implications of the non-trivial CP phase on the chargino 
sector. They also have computed the scalar-pseudoscalar mixings arising from the neutralino 
sector, and compared with those from the chargino sector [8]. Their calculations are quite 
exhaustive since almost all the relevant loops are included: top quark, the scalar top quark, 
bottom quark, the scalar bottom quark, tau lepton, the scalar tau lepton, W boson, the charged 
Higgs boson, the charginos, Z boson, the neutral Higgs bosons, and the neutralinos. They have 
reported that in the MSSM with explicit CP violation the neutralino exchange corrections to the 
mixings of the CP-even sector and the CP-odd sector are comparable to the chargino exchange 
corrections [8]. 

We would also investigate, in more detail, the Higgs sector of the MSSM with explicit CP 
violation at the one-loop level. In our investigation, we take explicitly the terms which are 
quartic in the electroweak coupling constants into account. We regard it consistent to consider 
the quartic terms in the one-loop effective potential contributed by the scalar top quark, the 
scalar bottom quark, the scalar tau lepton, the chargino, and the neutralinos, because the Higgs 
self contributions are essentially induced by them. 

Moreover, in our investigation, the pseudoscalar Higgs loop contribution is included, when 
the one-loop contributions to the neutral Higgs boson masses are evaluated, as well as all the 
loops of the relevant particles and superpartners. The inclusion of the pseudoscalar Higgs 
contribution is not only reasonable but also necessary, because there is a non-trivial correlation 
between the ordinary particles and the corresponding superpartners in the one-loop effective 
potential. In the neutral sector of the MSSM, there are four neutralinos as superpartners, and 
as ordinary particles there are Z boson and three neutral Higgs bosons, one of them being the 
pseudoscalar Higgs boson. Thus, the radiative contribution of the pseudoscalar Higgs loop should 
simultaneously be included, whenever the loops of Z boson and the two scalar Higgs bosons are 
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taken into account, corresponding to the loops of four neutralinos. But for the pseudoscalar 
Higgs contribution, the mass matrix for the neutral Higgs bosons at the one-loop level might 
be inconsistent and would yield incomplete masses for them and mixing angles among them. In 
order to obtain reliable results for the Higgs productions and their decays, it is quite necessary 
to include the pseudoscalar Higgs contribution. 

In this paper, we evaluate the masses of the neutral Higgs bosons in the MSSM with explicit 
CP violation at the one-loop level. We take into account the loops of the pseudoscalar Higgs 
boson as well as all the loops that have been considered in the investigations of Ref . [8] , namely, 
the loops of top quark, the scalar top quarks, bottom quark, the scalar bottom quarks, tau 
lepton, the scalar tau Icptons, W boson, the charged Higgs boson, the charginos, Z boson, the 
scalar Higgs bosons, and the neutralinos. We examine the effect of the contributions of the 
quartic terms and the pseudoscalar Higgs loops. We find that the mass matrix of the neutral 
Higgs bosons in the MSSM at the one-loop level is evidently affected by the inclusion of the 
quartic terms and the pseudoscalar Higgs contribution. 



2. The MSSM Higgs sector 

A starting point to evaluate the neutral Higgs boson masses in the MSSM with explicit CP 
violation may be the tree-level Higgs potential, which is given in terms of two Higgs doublets 
Hi (i = l,2) by 

y° = ^I{Hlmi+HlaH2r + ^-^{\H2f-m'f 

+ mllHif + ml\H2\'^ - ml{H'(eH2 + H.c.) , (1) 

where e is an antisymmetric 2x2 matrix with ei2 = 1, <? denotes the three Pauli matrices, gi 
and 32 are the U(l) and SU(2) gauge coupling constants, respectively, and {i = 1, 2, 3) are 
the soft SUSY breaking masses. Two soft SUSY breaking masses mf {i = 1, 2) may be assumed 
to be real, without loss of generality. We assume m\ = Imgje*'^'^ to be complex. 

Transforming to a unitary gauge, one may express two Higgs doublets in the mass eigenstates 

as 



Hi 
Ho 



1 I vi + hi + i sin f^hs \ 

J_( C0S/3C+ 

-^2 I V2 + h2 + icos Phs j ' 



where hi, h2, hs are three neutral Higgs fields, (7+ is the charged Higgs field, and is the 

relative phase (po between the two Higgs doublets. In the presence of the explicit CP violation, 
hi, h2, /i3 are mixed states of the CP parity. With respect to these three neutral Higgs fields, 
we have the corresponding minimum conditions as 

m\ = m|tan/3cos((;/!)3 -I- <^o) ^ cos 2/3, 

2 

m\ = m| cot /? cos(03 + (^q) H — ^ cos 2/3 , 
= m| sin(03 -|- ^o) , (3) 

respectively. 
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At the tree level, the phases (f)^ and (j)Q can be set to zero. Thus, it is well observed that 
there is no CP phase in the above tree-level Higgs potential. In order to accommodate the 
explicit CP violation scenario, the coefficients in the soft SUSY breaking terms are assumed to 
be complex in general. To be specific, we assume that the following quantities may be complex: 
The Higgs mixing parameters with mass dimension /i, the U(l) gaugino mass Afi, the SU(2) 
gaugino mass M2, and the trilinear soft SUSY breaking masses At, ^5, and A-,-. The complex 
phases in these quantities are responsible for the CP violation. On the other hand, we take the 
soft SUSY breaking masses mg, mr, mB, ttil, and be real. 

Now, as the electroweak symmetry is broken spontaneously, the neutral Higgs fields develop 
non-trivial vacuum expectation values (VEVs). The ratio of the two VEVs are defined as 
tan/3 = V2/V1. In terms of vi and V2, the fermion masses are given as rnf = [htV2)'^/2 for 
top quark, = {hbViY /2 for bottom quark, = [hrViY /2 and the gauge boson masses as 
f^w = (fl'2'y)^/4 and m| = [gl + gD^)"^ /'i with v"^ = vf + 

We calculate the radiative corrections to the tree- level Higgs sector. We employ the effective 
potential method [9] to estimate radiative corrections to the tree-level Higgs sector. The Higgs 
potential at the one-loop level is 

V = V^ + V^ , 

where represents the radiative corrections due to various loops, and is calculated by the 
effective potential method. It may conveniently be decomposed as 

= v^ + + + v^ + v'' + , (4) 

where 

3Mf / M? 
, 16^(i°Sa^ 

for the one-loop contributions of top quark and the scalar top quarks. 




327r2 I A2 2 167r2 \^ ^ A2 
for those of bottom quark and the scalar bottom quarks. 




A2 2 J 167r2 \^ A2 2 
for those of tau lepton and the scalar tau leptons. 





327r2 \^ A2 2 J 327r2 ^ A2 
for those of W boson, the charged Higgs boson and the charginos, 
, . 3M% /, Ml 3\ 4. Ml / Ml 3\ 

for those of Z boson and both scalar and pseudoscalar neutral Higgs bosons, and finally, 
4 TW^o / Mlo 

1=1 \ 
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for those of the neutralinos. 

In the above expressions, A is the renormahzation scale in the modified minimal subtraction 
(MS) scheme and Ai stands for the field-dependent mass matrix of the relevant ordinary particles 
and their superpartners. Note that the coefficient of differs from those of F* and as the 
color factor is absent in it. One may group the above components into two parts: V^, V^, , 
and into the charged part of , whereas F'' and into the neutral part. Note that the 
loops of all the relevant particles and superpartners are included: top quark, bottom quark, tau 
lepton, W boson, the charged Higgs boson, Z boson, the neutral scalar and pseudoscalar Higgs 
bosons, and then the scalar top quarks, the scalar bottom quarks, the scalar tau leptons, the 
charginos, and the neutralinos. 

In order to calculate the radiative corrections, one has to know the tree-level masses of 
particles and superpartners. The scalar top quarks, the scalar bottom quarks, the scalar tau 
leptons, and the charginos obtain their tree-level masses as follows: 



m 



tl,2 



rrit 

01,2 



m, 



t + \{mQ+ ml) + ^ cos 2/? ^ 



{rriQ - ml) 



So ^2 



cos 2/3 

,2 



+ (Aj + ix' cot^ /3 - 2At\i cot /3 cos ^t) 
ml + l{ml + ml) -"^008 2(3^ {^("^Q-^D 
^ cos 2/?! +ml {Al + jjk^ tan^ (3 



1 2 



2 ^biu tan /3 cos (j)},) 



ml + ^iml + mD-'^coslp^ {^(mi-m|) 



+ 



^2 ^2 

-mz - -mw 



m 



Xl,2 



^ cos 2/?| -I- ml{Al + ^? tw? (3 - 2At-h tan P cos (f>r) 
{^(M|-/x')-mif.cos2^|' 
-I- 2mw cos^ /?(M| + I? tan^ /? -|- 2M2/X tan /3 cos ^< 



(5) 



where D-terms are included. Notice that four non-trivial CP phases appear in the above ex- 
pressions. They arc defined as follows: (pt is the relative phase of At and )U, (f)^ is that of Ai, and 

ji, 0T- is that of A-T and /i, and (j)^ is that of M2 and ji. 
The tree-level neutralino mass matrix is given as 







M2 f,Hl 



\ 



91 ttO 92 rrO 



fP2 



-l^e^4>2 







(6) 



where additional complex phases appear: (pi is the relative phase between Mi and M2 and (f>2 
is the relative phase between M2 and /i. However, (/>2 is identical to cpc- Thus, the neutralino 
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mass matrix introduces only one more phase. Consequently, we have in general five non-trivial 
CP phases: (pt, (f'b, 4't, 4>c = <p2 and (pi. The above neutralino mass matrix is complex and 
symmetric, but not Hermitian. By diagonalizing the Hermitian matrix mIqM^o through a 
similarity transformation, the tree-level neutralino masses are calculated. They are denoted as 
m-o {i = 1, 2, 3, 4), and sorted such that m-q < "Z-o for i < j. 

For the Higgs sector, the tree-level mass of the pseudoscalar Higgs boson is obtained as 
2 _ 2ml cos(03 + ^o) 

Here we retain cp^ and (pQ because they become non-zero at the one loop level even though they 
are zero at the tree level. The tree-level masses of the charged Higgs boson and the remaning 
neutral Higgs bosons are given as 

2 2,2 



2 



1 



+ iTi^o =F Y (j^z + "^Ao) ~ ^m\rin?jifi cos^ 2(5 



(8) 



from the tree-level Higgs potential of the MSSM. Within the context of perturbation theory, we 
eliminate both and m| in the one-loop functions of the Higgs bosons by using the tree-level 
minimum equations. 

At the one loop level, there is one nontrivial mimimum condition with respect to /13 field. 
The CP-odd tadpole minimum equation is obtained as 

n 2 ■ / , , , X , 3r7T|/x^tshi0t 2 2 x , 3r7T|/xAsm^ 2 2n 

U = mo sm((i)3 + (Dq) -\ ;5 — tiimr , mr ) + ^ „ \ ^ tttt tn'mi 1 mi ) 

3 Y^ij i67r2'y2sin2/3-'^ *i' IQir^v^ cos^ P ^ ^'i' ^' 

ml^xAr s\n(pr 2 2 N m^l^M2 sin 0^ 2 2 n 
+ 167r2.2cos2/3 ^^^^"^^' W "^^^^ 



,_i47r2^;2 M A2 I nKo-m|o) ' 

'=-1 I \ / ; Xfc Xa 



where 



E3 = - (mto - M(){m%o - fi^)M2nm^ sin 02 

- (m|o - M|)(m|o - M^)Mi//(m| - mft.) sin(0i + ,^2) • (10) 

The dimensionless function fi is defined as 



i- / 2 2n 2 f 2, 2i o 

"^-^ ^ (mg - ml) A2 - - 2 ■ 

The six terms on the right-hand side of Eq. (9) come respectively from the tree-level Higgs 
potential, and the one-loop contributions of the scalar top quark, the scalar bottom quark, the 
scalar tau lepton, the chargino, and the neutralinos. 

By differentiating the Higgs potential at the one-loop level with respect to the three neutral 
Higgs fields, a 3 x 3 symmetric mass matrix Mjj is obtained in the (/ii, ^2, /i3)-basis. It may be 
decomposed as 
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where M^j is obtained from and M/^ from V^. In particular, M^^ may further be decomposed 



as 



where, as the superscripts suggest, Mfj represents the radiative contributions from F*. Likewise, 
Mfj, MJp mJ, Mj'j , and Mf- represent respectively the radiative contributions from V^, V^, 
V^, V^, and 

The matrix elements of Mij in the (/ii, /i2, ^3)-basis are easily calculated as 



Mil = rr?z cos^ 13 + fh\ sin^ [5 + Mii , 
M22 = m| sin^ P + m\ cos^ /3 + > 

mi + Mjg , 

— (m^ + m\) cos /3 sin /3 + M12 ; 



M33 
M12 

Mi3 

M23 



13 ) 



M 



(11) 



where the mass parameter uia is defined as 



m\ cos (^3 + ^0) + 



Sm^ /i A( cos (/); 



sin 2/3 

Sm^/iAfc cos (/!>6 
16lT^V^ cos^ /? 
m^/uM2 cos (/)( 



+ 



4 



47^2^,2 
2 

v 



-/i("i| 



IGTr^t;^ sin^ /3 

rn^fiAj- cos 



62- 



+ 



2 

XI' 



fe=l 



log 



A2 



IGtt'^v'^ cos^ /3 



33 



n (mto -mto) 



(12) 



with 



£^33 = (m^o - )(m?o - ix^)M2iim^ cos (/)2 

+ (m|o - Ml){mlo - lx'^)Min{ml - m^) cos(<^i + (^2) 



(13) 



In Eq. (11), the matrix elements Mj3 (z = 1, 2) represent the mixing between the scalar and the 
pseudoscalar components. Thus, Mj3 = M^g (i = 1, 2) indicate that the mixing occurs only at 
the one loop level. There is no mixing at the tree level. 

For M^., we first calculate M/J- from V^. The result is given as follows: 



Mi\ 



nir. 



+ 



+ 



327r2 
m% cos^ /3 



[m 



log 



^ H 

mf^ cos' j cos- 2/3 , 

^^22 ^log 

32yr^t;^ 

,2 «a2 



,2\2 




327r2?;2 



(m^o sin^ /? — 4m| cos^ /3)/i(m^, m|^) 



m| cos^ l3Ahi log{m1/m\) 



167r2 



(m|^- 



m'y cos^ /3 / m% 
' log -r# 



87r2-t;2 



M 



22 



sin^ /3A^^ /2(m2, m\) m| cos^ /3 



327r^ {rnjj — m 



2\2 



327r2t;2 



m 



AO 



^ A6 ^ 

4m|)/i(mL m%) 



7 



M. 



33 



M 



12 



m\ SVC? (3 



+ 



+ 



8n^v 





2sm^2/3-3)}/iK> "^?/) 



m^m^o cos^ 2/3 




647r2 {m\j — mfj'^ 



+ 



m 



|sm2^ 2 



327r2i;2 



-(m 



AO 



171% sin 2/3 



log 



m| sin 2/3 cos^ 2/3 



+ 



• sin 2/3 



647r2 



(A;,, + A;,, 



g4y^2^2 



log 



A" 



A2 



+ 



m. 



■ sin 2/3 



-2m|)/i(mL ml,) 

log (mf,/m|) 
(m|, - ml) 



rrir 



Mfa = 0, 



M. 



23 



(14) 



where 



/2(m^, m ) = — ^ log ^ 



m 



y 



mt 



mt 



A., 



2m 



2 \ 2 , 4m^m^o . 2 o 
m.ojcos p -\ sm /3 , 



2mi 



{mz - m^o)sin /3 + 



'AO 



cos^ /3 . 



(15) 



In the above expressions, m/jO, and rrifjo, and m^o are the tree-level masses of the scalar and 
the pseudoscalar Higgs bosons. Since there is no CP phase in V'^, we have = iVf^s = in 
Eq. (14). Therefore, the explicit CP violation scenario in the Higgs sector of the MSSM may be 
regarded as the radiative CP violation, in the sense that the CP violation occurs only through 
the radiative corrections due to superpartners, via the CP phases in soft SUSY breaking terms. 

Next, we calculate M<^^ , the contributions of radiative corrections due to the four neutralinos. 
We obtain 




A2 I dhi 




where the first-order derivative dmio/dhi is given explicitly by 



dm?-o 
dhi 



Aimlo + BiVnU + QmL + A 

A/j. Afc A/j 



n(m|o 
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and the second-order derivative d^niio / dhidhj by 

x" _ _ ■' xj ■' XI ■'XI ■' 
a^k Xo. 

I I -^fc Xa I Xa Afc I n?"! 

^^^^(m|o-m|o) dh, dhj ^ dhi dhj j ■ ^ > 

The formulas for coefficients Ai, Bi, Ci, and Di (i = 1, 2, 3) are given in Appendix A, and those 
for Aij, Bij, Cij, and Dij, = 1,2,3) are given in Appendix B. Actually, Aij = for all i,j 
because four elements of the neutralino mass matrix are zero. 

We also carry out calculations for the remaining components of the one-loop mass matrices, 
namely, M/^, Mfj, M[j and M^j, and the results are given in Appendices C, D, E, and F, 
respectively. 

In the scenario of explicit CP violation, among the elements of the mass matrix for the 
neutral Higgs bosons, Mis = M}^ = M^^ + Mfg + + M^l + M^l {i = 1, 2) are responsible for 
the scalar-pseudoscalar mixing and eventually account for the CP violation. It is worthwhile to 
notice that our calculations show the general dependence of Mjs on the various CP phases: M/g 
{i = 1,2) are proportional to sin^j, M^^g (z = 1, 2) to sin (i = 1, 2) to sin0^, {i = 1,2) 

to siiKpc. Meanwhile, {i = 1,2) arc zero, as noticed above. As for [i = 1,2), they have 
complicated expressions, but one can easily sec that every term is proportional to sin(0i + <p2) 
or sini?^! since B^, C3, D^, are proportional to them whereas B-i-^, Ciz, D^^ (i = 1,2,3) are zero. 
Consequently, if any of the CP phases be not zero, they might give rise to the CP violation via 
the scalar-peudoscalar mixing at the one-loop level. 



3. Numerical analysis 

We investigate the effects of CP phases in the one-loop corrected MSSM Higgs sector. As we 
have emphasized, if there exist some CP phases in the Higgs potential at the one-loop level, they 
would produce the scalar-pseudoscalar mixing through non-zero M13 and M23 and thus generate 
explicit CP violation in the neutral Higgs sector. In the MSSM Higgs sector, we show that there 
are in general five independent CP phases. They are contained in the masses of the scalar top 
quarks, the scalar bottom quarks, the scalar tau leptons, the charginos, and the neutralinos. 

The mass matrix for the neutral Higgs bosons is derived by differentiating twice the Higgs 
potential with respect to the three Higgs fields. It is expressed as a 3 x 3 matrix. By calculating 
its three eigenvalues, one can obtain the three masses of the neutral Higgs bosons, analytically 
in principle. However, we would not write down the full expressions for the neutral Higgs boson 
masses, since they are very complicated functions of many parameters coming from the soft 
SUSY breaking in the MSSM. Rather, we would like to analyze numerically the effect of the CP 
phases on the mass matrix of the neutral Higgs boson by focusing our attention to M^a {i = 1,2). 

For the numerical calculations, we need concrete numbers for the parameter values. We set 
mt = 175.0 GeV, = 4.0 GeV, = 1.7 GeV, mw = 80.4 GeV, mz = 91.1 GeV, and sin^ Ow 
= 0.231. Then the remaining relevant free parameters are A, tan/3, fj,, mq, uit, mB, ttil, mE, 
At, A},, At, Ml, M2, 4>t, 4>b} (pT, 4>c{4>2), and 4>i. Since for a moderate values of the tan/3, the 
contribution for the scalar bottom sector as well as the scalar tau lepton one is relatively small, 
we set niQ = rriL, mx = ruB = rriE, At = Ai, = A^-. At the electroweak scale one can take 
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the relation of Mi = 5tan^0^yM2/3 between U(l) and SU(2) gaugino masses. Thus, we have 
A, tan/3, tuq, itit, At, M2, and five CP phases as free parameters. Note that we will choose 
different values for the soft SUSY breaking masses of SU(2) doublet and singlet scalar fermions 
such that niQ = is different from uit = tub = tue in order to consider large mass splitting 
between the left- and right- handed scalar fermions. 

Now, as a typical illustration, we set <pt = 4>b = (pr = 4'ci4'2) = 01 = vr/3. We fix the renor- 
malization scale A in the effective potential at 300 GeV. We set the remaning free parameters 
as fhA = 300 GeV, mq = 800 GeV, mr = M2 = 400 GeV, and At = 200 GeV. With these 
input values, we calculate the elements of the neutral Higgs boson mass matrix for some values 
of tan (3 and /i. 

Our results are shown in Tables 1-8. These Tables are divided into two categories: in Tables 
1-4 the contributions from the quartic terms of 0{gf) {i = 1, 2) and the pseudoscalar Higgs 
contribution are excluded whereas they are included in Tables 5-8, in order to exhibit the effect 
of these contributions. Tables 1, 2, 5 and 6 displays numerical results for small tan (3 = 5 while 
Tables 3, 4, 7 and 8 for relatively large tan/3 = 30. For Tables 1, 3, 5 and 7 we set = —400 GeV 
while for Tables 2, 4, 6 and 8 we set /it = 400 GeV. The entries of Tables are the elements of the 
mass matrix of the neutral Higgs bosons in the {hi, h2, /i3)-basis, where all the contributions are 
fully listed. The unit is (GeV)^. Note that, as the CP phases induce the CP violation through 
the scalar-pseudoscalar mixing, the columns for M13 and M23 are filled with generally non-zero 
numbers. The number in the last row in each column is the sum of the numbers in the preceding 
rows in the same column. 

In all of the Tables 1-8, one can easily notice that, among the contributions to M13 and 
M23, the tree-level ones and the Higgs ones vanish: M^^ = M^^ = and Mfg = M^^ = 0. 
This reflects the fact that there is no CP phase in the tree-level Higgs sector of the MSSM. 
All the other elements are not zero, accounting for the scalar-pseudoscalar mixing, because the 
CP phases are present there. Actually, one can trace the effects of the CP phases on various 
contributions to M13 and M23. If (pt = 0, one would have M^^ = M23 = 0. Further, one would 
have = M23 = for (pi, = 0. Likewise, one would have Mi^ = M23 = for (pT = 0, 
M^g = = for (pci(p2) = 0, and M^g = = for (j)2 = (pi = 0. Among the non-zero 

contributions to M13 and M23, we find that M^^ of Table with = —400 GeV contribute about 
25 ~ 30 % to the scalar-pseudoscalar mixing. Thus, Tables 1, 3, 5, and 7 indicate that the 
radiative corrections due to the neutralino loops contribute roughly 25 % to the CP mixing 
between hi and /13 components. 

Now, let us compare Table 1 with Table 5, Table 2 with Table 6, and so on to examine the 
effects of the contributions from the quartic terms of 0{gf ) {i = 1,2) and the contribution by 
the pseudoscalar Higgs boson A^. It is quite evident that there are non- negligible differences 
between them. The effects can most clearly be seen in M33. In Tables 1-4, we see that M33 is 
identically zero. This is expected since neither A^ nor quartic term contributions are included in 
Tables 1-4. On the other hand. Tables 5-8 show M33 is definitely far from zero, indicating that 
both of them contribute in Tables 5-8. For the particular parameter values, -M33 is calculated 
as 55.9, 58.0, 30.2, and 39.2 (GeV)^ in Tables 5,6,7, and 8, respectively. 

The effects of the contributions from the quartic terms and A^ are also seen in other entries, 
too. By comparing Table 1 with Table 5, for example, one can see that those contributions 
alter the values of M{\, M^^^ and M^^- The value of M3'*3, as well as those of M^\, M^2^ and 
Mi2- is essential to evaluate the neutral Higgs boson masses from their mass matrix and the 
mixing among them, at the one-loop level. Consequently, these contributions should be included. 
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as they play roles definitely in the matrix elements M-^j at the one-loop level. Without these 
contributions, the neutral Higgs boson masses at the one-loop level would be changed, if the 
changes be small. We further note that M33 is in particular also affected by the inclusion of 
these contributions. 

Now, to be concrete, we evaluate the neutral Higgs boson masses. From Table 5, we obtain 
mh,{i = 1,2,3) = 109.8, 300.1, 302.6 GeV and we obtain 108.6, 300.1, 302.8 GeV from Table 
6. Since Tables 5 and 6 differ in the value of n, the dependence on jj, is relatively insignificant, 
for tan/3 = 5. For relatively large tan/? = 30, the masses of the three neutral Higgs bosons are 
obtained as mh^i = 1,2,3) = 115.2, 300.0, 300.1 GeV for Table 7 = -400 GeV) and 115.0, 
300.0, 300.2 GeV for Table 8 = 400 GeV). We see that the mass of the lightest neutral Higgs 
boson becomes slightly larger for large tan /?. All of Tables 5-8 produce the neutral Higgs boson 
masses without contradicting LEP2 data. 

We have some comments on the radiative corrections due to the loops of Z boson, the neutral 
Higgs bosons, and the neutralinos. One can observe in Table 5-8 that the sum of AI12 and Mj^g 
are comparatively large. This implies that the contributions of Z boson, the neutral Higgs 
bosons (both scalar and pseudoscalar) , and neutralinos play important role in the radiative 
corrections to the mixing between hi and /12 components at the one-loop level. Especially, this 
mixmg occurs most largely in Table 6 as {M^2 + ^1*2 ) = -78.1 (GeV)^. 

The contributions of the neutralino loops depend crucially on the CP phase (pi. In other 
words, the CP phase (pi occurs only in the expressions for the neutralino contributions. Now, 
in order to examine in more detail the dependence of the contributions of the neutralino loops 
on (pi, we plot , M13, M23 , and M23 as functions of (pi in Fig. 1, where the values of other 
parameters are the same as Table 5. We do not plot M/3, M^^ and M^^ and M^^ in Fig. 1, 
since their contributions are the same as those numerical values of Table 5. In Table 5, the 
contributions of the loops of the scalar quark and the scalar tau lepton are almost null for M13. 
For Mi3, the contributions of neutralino loops are larger than those of the scalar fermions for 
some values of (pi . For both M13 and AI23 , Fig. 1 shows that the contributions of the neutralino 
loops are smaller than the chargino ones for the whole range of < 0i < 27r. 

4. Conclusions 

In the MSSM, the tree-level neutral Higgs sector may be divided into the scalar part and pseu- 
doscalar part, and there is no mixing between them. Any phase that can cause the scalar- 
pseudoscalar mixing, hence the CP violation, can be absorbed away at the tree level. At the 
one-loop level, where explicit CP violation is viable by introducing several CP phases in the effec- 
tive Higgs potential, the scalar-pseudoscalar mixing occurs in general. The scalar-pseudoscalar 
mixing is manifested by the non-vanishing M13 and M23 matrix elements of the neutral Higgs 
bosons. Evidently, these off-diagonal elements affect the masses of the neutral Higgs bosons 
when the mass matrix is diagonalized. 

The mass matrix of the neutral Higgs boson in the MSSM is evaluated at the one-loop 
level with explicit CP violation. In explicit CP violation scenario, five non-trivial CP phases 
are introduced, from the soft SUSY breaking terms of the MSSM Lagrangian, in the masses of 
the scalar top quarks, the scalar bottom quarks, the scalar tau leptons, the charginos, and the 
neutralinos. These phases penetrate into the mass matrix of the neutral Higgs bosons. All the 
contributions of relevant loops are taken into account: The loops of the pseudoscalar Higgs boson 
as well as all the loops of top quark, the scalar top quarks, bottom quark, the scalar bottom 
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quarks, tau lepton, the scalar tau Icptons, W boson, the charged Higgs boson, the charginos, Z 
boson, the scalar Higgs bosons, and the ncutralinos. 

Especially, we consider the contributions of the terms quartic in the electroweak coupling 
and the pseudoscalar Higgs loop contribution in the one-loop effective potential to the neutral 
Higgs boson masses at the MSSM with five explicit CP phases. It is found that the contributions 
of quartic terms and are definitely non-zero to the (3,3)-element of the mass matrix of the 
neutral Higgs bosons. The CP mixing between hi and hs components can be induced about 25 
% of the considered total contribution by the neutralino contributions. The contributions of Z 
boson, the neutral Higgs bosons (both scalar and pseudoscalar), and the neutralinos contribute 
largely to the mixing between hi and /12 components above other ones. 
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Appendix A 

The coefficients that appear in the first derivatives of the neutralino masses with respect to the 
neutral Higgs fields in the radiatively corrected mass matrix for the neutral Higgs bosons are 
given as follows: 

4m 7 cos d 
Ai = ^ , 

V 

AMlmL cos f3 4M| (m| - mL) cos (5 4M| (m| + ji^) cos (3 
Bi = 1 1 

V V V 

4M2in^fi sin /? cos ^2 4:Mifi{m^r — rn^) sin /? cos(<^i + (^2) 
I ) 

V V 

4m|//^ sin/3sin2/3 4Mf mfj/(m^ + cos /? 

V V 

^ AM^{Ml - ml,){ml, - m| - /x^) cos P 

V 

^ SMxM2m'^{m^ - m|)cos/3cos^i ^ 4:M2m^ii{Ml + /x^) sin/? cos 02 

V V 

4.Miii{m% - 'mw){M^ + M^) sin/?cos(0i + ^2) 

V 

AMlm^lx^ sin/?sin2/3 4.M^iJ?'{m% - f^wf sin/3sin2/3 
Di = H 1 

V V 

^ 8MiM2m^/U^(m| - m^^) sin /3 sin 2(3 cos 0i 4Mf M2m^/Li^ sin /3 cos ^2 

f V 

+ , (18) 



and 



and 



A2 = ^1 (cos /3 •i-^ sin /3) , 

B2 = -Bi(cos/9 sin/?) , 

C2 = Ci (cos /3 sin /?) , 

1)2 = £'i(cos/3^sin/3) . (19) 



^3 = 0, ^ 

4M2)»^7v sin (:;2 4i\/i/((»;^' — irr^ ') H'm{pi +02) 



B 



3 



^ 4M2m^//(Mf + //2)sin(/)2 , 4Mi//(to| - to2^)(M| + gj^^^^^ ^ ^2) 
1-^3 = r 



■u V 
^ 4Mf Mamf^M^ sin (/)2 ^ 4Mi M2V(mft. - m|) sin((/)i + (^2) ^20) 



Appendix B 

The coefficients that appear in the second derivatives of the neutralino masses with respect to 
the neutral Higgs fields in the radiatively corrected mass matrix for the neutral Higgs bosons 
are given as follows: 

An = 0, 
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and 



and 



Bl2 
Cl2 



8m% cos^ /9 

SM^m^ cos^ p 8M|(m| - mfj,)^ cos^ /? 
16MiM2m|(^(m|^ — m|) cos^ /3cos 0i 

~^ 2 ' 

L>n = 0, (21) 



Bii 
Cii 



^22 = ^11 (cos /3 sin /?) = , 

-B22 = -Bii(cos/3 sin/3) , 

C22 = Cii(cos/3 ^ sin/3) , 

D22 = Du{cos 13^ sin/3) = , (22) 



A12 = , 

4m^ sin 2/3 



^2 



4M^m^ sin 2/3 8m|/x2 sin 2/3 4M| (m| - m^)^ sin 2/3 

^2 y2 y2 

8MiM2m^(m^ - m| ) sin 2/3 cos ^1 
2 ' 

SMfm^n'^ sin 2^ SMf/x^ (m| - m^)^ sin 2/3 

-L'12 = n 1 5 

16MiM2mft^//^(m| - m^) sin 2/3 cos (/>! 

+ ^;2 ' '^^"^^ 

and = Bi3 = = A3 = for Z = 1, 2, 3. 



Appendix C 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the top quark and scalar top quarks are 

_ 3 [ m,VA,-^ cos/3A,- | V2(m|, m|) Smj cos^ /? / ^ ^| ^^ 

-"^^11 - ~ ^_2„.2 1 „;„ a "I n r 7ri2 Z:2~;2 ^ c^_2„.2 
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+ 



3 cos^ /3 / 4m^ 




Att^v'^ I sin/3 2 J (mj -mj y GAir'^v'^ \ 

2 

w 



^ 3m| cos/3 irnh^ ^ cos/3A,-\ log(m^-ym^-^ ^ 
Btt^v^ I sin/3 



3 m|^tA 




47r2 



2yltA,-^ sin/3A,-\ log(m4/m|J 



Btt^u^ I sin^ /3 / I sin /3 2 ( (m| — m? ' 
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^33 



3 / 2m| sin /3m- 
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with 




3m| sin 2/3 
1287r2?;2 \^sin2/3 

47r2?;2sin2/3 \ sin/3 ' 2 
_^ Sm'^mlfiAt cot /3 sin log(m|/m| ) 
167r27;2 sin /3 (m? — m? ) 



|. 

23 - 47r2?;2 sin2 /3 1 sin/3 2 / (m? - m2 ^2 



, _ Sm^tl^At sin \ ^iM^i^ sin /3A,- 1 h {ml , mf^ 



t2 tl^ 

SmIfiAt Sin 4>t ( Jm^ _ 2 \ ^o^Hjml) 
^ 167r2t;2sin/3 \s\n^ P j {m]^-ml) ' ^ ' 



A^^ = /X cot (5 — At cos (?!>t , 



= At — II cot /3 cos ( 

Af = Q'Ti^ - ^'Til^ - mf. + - ^"t-I) cos 2/3| . (25) 



Appendix D 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the bottom quark and scalar bottom quarks are 

M'u = - 
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^23 



with 




-^22 - - -JZ^ 1 -3 ( TZ2 -^^^ + 54^2^2 1 ^4 



^ 3m| Bin /? f mgy Aj^ sin /j A, ] logK, /"j, ) 
87^2^2 I 2 I (m? — m? ) 

h _ omi^fi .i^sm 9fj''-v /^^ 



02 01 ' 




IQtt'^v'^ \ cos^ /3 / I COS 6 2 f (m? — m? ) 

\ /V. ^ 02 01 ' 

327r2i;2 

3m| sin 2/3 
1287r2t;2 

Smf^^ft sin 
47r2'u2 cos2 /3 

_^ 3m2^yl6sin(/)5 



+ 



cos /? 

J _ 3ml fi Ah sin 4>b 




Att'^iP cos'^ 3 I cos/3 2 ( (m? — m? )2 

V V ^ 62 Oi' 

^ 3m2m|^Abtan/3sin</)fclogK,/^t 

IGtt'^v'^ cos /3 (m? — TO? ) ' 

^02 01 ' 



Ag^ = — /Lt tan /3 cos , 
Afoj ~ ^ ^™ P — Ah cos </>b , 

Aft = - I'^'w^ {"^Q - + Q"i| - ^"^w) cos 2^1 . (27) 



Appendix E 



The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
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of the tau lepton and scalar tau leptons are 



^ 1 j mlArAf^ ^ cos/3A^ 



M. 



with 




/If" 

47r2-t;2 cos4 /? (m?^ - 

_ 1 [ m^A^Af, ^ cos/3A^ | j rn^fiAf^ smpAf \ f2{m\, m%) 



sin/3Afl log(r7T,| /m? 




(m4 -m|J 



Af ] l og(m4/^fi) 
(m? - m? ) 



cos/3Afl /2(m? m?) 



13 - 47r2^;2cos2/? \ cos/3 2 j (m?^-m?J2 

^ ml fiAr sin (pr f 2 ^ log ("^12/"^^ 



^ _ mliiAr sin (br { ml^ A sin /3 Af 1 /2 , rn? J 



23 4^2^2pQg2^ [ cos/3 2 ] {ml^-m\, 



m^m^iiA^ tan /3 sin log (m?^/m? J 
IGtt^v^ cos /3 (ml^ — ) 



Af^ = A-T — 11 tan /3 cos 0^- , 
Afj = /X tan P — At cos , 
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= - -ml+ Qm| - m^^ cos 2/3| . (29) 

Appendix F 

The elements for the mass matrix of the neutral Higgs bosons due to the radiative contributions 
of the W boson, charged Higgs boson, and charginos are 

a4 = g^(4„;,A^. A,. - A,)^ f <r^- ' f ;^ + log ( 



with 



XI X2 



cos^ /? (4m^M2A^, - A^) / \ cos^ [3 ^ 



-/i("^xi 



2 2 



-^22 = 4mfv.M cot /?A^2 + A^ % 2 TT? + Too log 4 4 

^ X2 Xl' \ Xl X2 / 

sin^ /3 (4m^/x cot pA^, + A^^) /m|^ \ m^sin^/j o ^2 ^ 

Zri- o_o_.o /U"^xi. "^yJ ' 



2(;»2,-/;il/2siiiaJ-,/2('"|,- '4J 



X2 XI 

^12 = t^^(4M2.M2A^, -A^)(4m2^/v,cot/3A^2+A 



sin2/3,,,^2 A A N/.„„„2 «A , A ^/2(rn|^, 



16^2^2 V*--W/--2--xi --XA-'-iyP --'^ A^--X2 -r --X7 ^^2 _ ^2 )2 

cos P (M2 sin /3A^^ + ft cos /JA^^) ^ / m^^^ \ 

TT^V^ ("^X2-"ix2) °^V"^xJ 



m^sin2/3 m\yml,+ \ m^sin2/3 2 2 

+ log (^7^5;^; J + ^;^/i(-x.> -X 

MX _ m^M2MCOs/3sin<^e (4m^M2A^, - A^^) , ^2 

m^M2;U cos /3 sin 
7r2t;2(m|2 — 

^ m^M2/xsin/?sin(/.c (Amir cot PA^^ + A^) ^ 2^ 

,4 




my^M2pi sin /? sin <;i!), 
Tr^v^{m^^ - m|J 




(30) 



Aj^^ = M2 + /U tan P cos 
A^2 = -^2 cos 4>c + IJ- tai 



M2 COS 4>c + IJ- tan /? , 

Aj^ = 2mvi/(M| - /x^ - 2mvi.cos2/3) . (31) 
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Figure Captions 

Fig. 1 : The plot of the (1,3)- and (2, 3)-elements of the mass matrix of the neutral Higgs 
bosons at the one-loop level with the contribution of terms of 0{gf) {i = 1, 2), M^^ (dot-dashed 
curve), (dashed curve), M13 (dotted curve) and M23 (solid curve), as a function of (f>i for 

= = = 0c(<^2) = 7r/3, A = 300 GeV, and rfiA = 300 GeV,mQ = 800 GeV, mx = 400 
GeV, At = 200 GeV, and M2 = 400 GeV. We set tan /3 = 5 and /x = -400 GeV. These values 
for the parameters are the same as Table 5, except that is taken as a variable. 
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Table Caption 



Table 1: The elements of the symmetric mass matrix of the neutral Higgs bosons in the 
{hi, /i2, /i3)-basis, at the one-loop level with CP violation for (pt = (pb = (pr = 4>c{4>2) = <pi = 7r/3. 
Here, the contributions of terms of 0{gf) {i = 1,2) and are neglected. The unit is (GeV)^. 
The values of the relevant parameters are A = 300 GeV, ffiA = 300 GcV, mq = 800 GeV, rriT 
= 400 GeV, At = 200 GeV, and M2 = 400 GeV. We set tan/3 = 5 and = -400 GeV. The 
number in the first row in each column is the tree-level value. The number in the last row 
in each column is the sum of all numbers in the preceding rows, representing the value at the 
one-loop level. The numbers in each column in between the two rows represent the various loop 
contributions as decomposed in the one-loop Higgs potential. 







(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


n ,r() 

«j 


86856.9 


11424.0 


90000.0 


-18900.1 


0.0 


0.0 


IVl ■ ■ 
11 


— 7.2 


6362.6 


-6.9 


229.3 


7.0 


-225.0 


iVl ■■ 
11 


0.1 


-0.008 


-0.001 


0.01 


-0.007 


0.003 




0.002 


-0.0001 


-0.00001 


0.0002 


-0.00009 


0.00004 




3.0 


-268.3 


16.5 


-31.5 


17.1 


40.7 




-3.7 


-150.4 


0.0 


-12.3 


0.0 


0.0 




0.1 


-146.1 


12.4 


-17.6 


10.2 


31.3 




86849.3 


17221.8 


90022.0 


-18732.3 


34.4 


-152.9 




Table 2: 


The same 


as Table 1, except for fi 


= 400 GeV. 




(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


M". 

11 


86856.9 


11424.0 


90000.0 


-18900.1 


0.0 


0.0 


11 


-0.08 


6267.7 


-6.9 


-25.5 


0.7 


226.0 


M^- 

'11 


0.1 


-0.006 


-0.001 


0.01 


0.009 


-0.003 


ij 


0.002 


-0.00008 


-0.00001 


0.0002 


0.0001 




-0.00003 


M^. 

11 


-14.6 


-288.7 


16.8 


-77.9 


1.7 


-37.1 


M^- 

11 


-3.6 


-147.2 


0.0 


-13.2 


0.0 


0.0 




-9.3 


-155.1 


12.7 


-42.1 


-1.9 


-29.9 




86829.3 


17100.5 


90022.6 


-19059.0 


0.6 


158.9 
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Table 3: The same as Table 1, except for tan (3 = 30. 



{^.3) 


(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


M% 


89909.3 


8371.6 


90000.0 


-3272.3 


0.0 


0.0 




-2.7 


6088.2 


-6.4 


139.1 


4.1 


-212.8 




7.2 


-9.3 


-1.7 


1.7 


-0.7 


4.0 




0.1 


-0.1 


-0.02 


0.02 


-0.009 


0.05 




6.7 


-279.2 


16.7 


13.7 


10.9 


40.0 




-0.1 


-170.2 


0.0 


-1.6 


0.0 


0.0 




2.1 


-155.7 


12.5 


7.1 


5.3 


31.3 




89922.7 


13845.2 


90021.1 


-3112.2 


19.7 


-137.3 


Table 4: The same as Table 2, except for tan /? = 30. 




(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


<• 


89909.3 


8371.6 


90000.0 


-3272.3 


0.0 


0.0 




-1.6 


6073.0 


-6.4 


-106.5 


3.2 


212.9 




3.2 


-9.0 


-1.7 


6.2 


2.7 


-3.9 




0.05 


-0.1 


-0.02 


0.08 


0.03 


-0.05 


t7 


3.7 


-309.5 


16.7 


-33.0 


8.3 


-39.4 


-'"J7 


-0.1 


-136.4 


0.0 


-3.0 


0.0 


0.0 


1^ 

17 


0.5 


-157.3 


12.6 


-17.4 


3.2 


-31.1 




89915.1 


13832.2 


90021.2 


-3426.1 


17.5 


138.4 
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Table 5: The elements of the symmetric mass matrix of the neutral Higgs bosons in the 
(/ii, /i2, ^3)-basis, at the one-loop level with the contributions of terms of 0{gf) [i = 1,2) and 
A^. The setting of the parameters is the same as Table 1. 





(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


u 


86856.9 


11424.0 


90000.0 


-18900.1 


0.0 


0.0 




-0.7 


5832.3 


-7.0 


265.7 


4.3 


-212.3 




1.0 


30.6 


-0.001 


-5.1 


0.02 


-0.1 


MT. 


0.1 


4.7 


-0.00001 


-0.8 


0.002 


-0.01 


y 


1.7 


-311.8 


16.6 


-23.8 


17.3 


40.0 




53.1 


-148.6 


55.9 


-31.5 


0.0 


0.0 


17 


2.8 


-82.3 


12.4 


-15.5 


10.5 


31.8 




86915.1 


16749.0 


90069.7 


-18711.4 


32.2 


-140.6 


Table 6: The same as Table 5, except for ^ = 400 GeV. 




(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 




86856.9 


11424.0 


90000.0 


-18900.1 


0.0 


0.0 




0.1 


5744.0 


-6.9 


24.7 


3.5 


213.3 


11 


1.0 


30.5 


-0.001 


-5.3 


-0.02 


0.1 


W- 

11 


0.1 


4.7 


-0.00001 


-0.9 


-0.002 


0.01 


y 


-16.2 


-331.7 


16.9 


-69.6 


1.6 


-36.4 


I- 3 


53.9 


-146.1 


58.0 


-31.8 


0.0 


0.0 


4^ 

11 


-6.7 


-95.5 


13.0 


-46.3 


0.5 


-30.1 


M,j 


8C889.3 


1GG30.0 


90072.5 


-19029.5 


5.G 


116.9 
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Table 7: The same as Table 5, except for tan (3 = 30. 





(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 




89909.3 


8371.6 


90000.0 


-3272.3 


0.0 


0.0 


M^- 
11 


-2.1 


5561.2 


-6.4 


139.2 


3.7 


-200.4 


It 


6.7 


51.3 


-1.7 


8.4 


-0.5 


-2.4 




0.07 


7.0 


-0.02 


0.3 


0.006 


-0.4 


'ti 


6.7 


-324.4 


16.7 


14.7 


10.9 


39.2 


'11 


40.0 


-168.3 


30.2 


-4.0 


0.0 


0.0 




5.0 


-91.3 


12.8 


12.8 


8.1 


32.1 


Mi,- 


89965.7 


13407.2 


90043.2 


-3100.9 


22.4 


-131.8 


Table 8: The same as Table 6, except for tan /3 = 30. 




(1, 1) 


(2, 2) 


(3, 3) 


(1, 2) 


(1, 3) 


(2, 3) 


<• 


89909.3 


8371.6 


90000.0 


-3272.3 


0.0 


0.0 


M*. 
11 


-2.0 


5547.0 


-6.4 


-92.1 


3.6 


200.5 




3.2 


51.1 


-1.7 


5.5 


2.4 


2.4 




0.06 


6.9 


-0.02 


-0.1 


0.01 


0.4 


ti 


3.6 


-354.7 


16.8 


-31.2 


8.3 


-38.6 




38.9 


-140.8 


39.2 


-3.8 


0.0 


0.0 


11 


2.5 


-93.6 


12.9 


-22.7 


6.0 


-31.7 


Mij 


89955.7 


13387.7 


90052.3 


-3417.0 


20.5 


133.0 
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Fig. 1: The plot of the (1, 3)- and (2, 3)-eIements of the mass matrix of the neutral Higgs bosons 
at the one-loop level with the contribution of the quartic terms of 0{gf) (i = 1,2), M^^ (dot- 
dashed curve), M^^ (dashed curve), Mi3 (dotted curve) and M23 (solid curve), as a function of 
(pi for (j)t = = (j)^ = (t)c{(t)2) = vr/3, A = 300 GeV, and niA = 300 GeV, mq = 800 GeV, mr 
= 400 GeV, At = 200 GeV, and M2 = 400 GeV. We set tan/3 = 5 and /x = -400 GeV. These 
values for the parameters are the same as Table 5, except that 0i is taken as a variable. 
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